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ABSTRACT 


This paper presents a practical statistical overview of the pressure- 

| volume calibration curve for large nuclear materials processing tanks. 
It explains the appropriateness a: SRR splines (piecewise 
polynomials) to this curve, and it presents an overview of the 
associated statistical uncertainties. In order to implement these 
procedures a practical and portable FORTRAN IV program is provided along 
with its users' manual. Finally, the recommended procedure is 


demonstrated on actual tank data collected by NBS. 


Key Words: Volume calibration; differential pressure; splines; 
accountability; statistics. 
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1. Introduction. 

The direct measurement of liquid volume in large processing tanks, 
especially with internal structure, is impractical at best. Measuring 
(differential) pressure is simple and quick. However, in order to 
estimate the volume indirectly by observing pressure, it is necessary to 
use the relationship between volume and pressure. This relationship is 
known .as..a calibration curve; its estimation is the process known as 
calibration. 

Fitting a calibration curve is much like regression, in that for 
"known" values vy of volume, one obtains one or more observations Pij of 
the. differential pressure py = p(vy), and "fits" a response function 
p(v) .by statistical methods - usually by least squares. At this point, 
the correspondence stops. Whereas regression is used to predict values 
of the dependent variable (p) for given values of the independent vari- 
able (v)...or to. test a proposed relationship between the variables, a 
calibration curve is used to estimate values of the independent variable 
v corresponding to new measured values of the dependent variable 
pe Furthermore, the confidence interval (or uncertainty measure) which 
is desired is not for p, but rather for v. And finally, systematic 
error ie intxaduced by tae of fit of the calibration curve, and ina 
materials eahdnn coud situation this may be crucial. 

THRE goceuanEesinte method for producing valid uncertainty limits 
for Beet creole badkued libracton curve by using calibration 
Eeueuions eiiichoatemsnaoths piecewise polynomial functions called 
“spldmess" Taking advantage of an ae ales to calibration originated by 


Scheffé [1] and further elucidated by Scheffé, Rosenblatt and 


Spiegelman [2], it provides statistically sound uncertainty limits, not 
just for a single estimated value of volume, but for all volumes 
estimated by use of the fitted curve. This approach overcomes a major 
theoretical problem with earlier methods: it makes proper allowance for 
the contribution to the overall uncertainty of errors in fitting the 
curve. 

The procedure presented herein has been implemented, based upon a 
spline-fitting program due to deBoor [3]. The resulting FORTRAN program 
has been tested on various sets of data, including actual tank data. 

The remainder of this paper is organized as follows. Section 2 
contains a discussion of the pressure-volume model, and the statistics of 
calibration. Section 3 contains a discussion of an example, and of the 
printout produced by the program. Section 4 is essentially a users' 
manual for the program. . In Section 5 will be found a discussion of 
open questions, work in progress, cautions, and possible extensions of 
this technique. Finally, Appendix 1 contains a discussion of Scheffe's 
constant c, used as an input to the program, and Appendix 2 contains a 


listing of the program. 


2. The Scientific Basis for, and Interpretation of, a Calibration Curve. 

In order to provide statistically valid uncertainty limits for the 
volume estimates obtained through the use of a calibration curve, it is 
necessary to have a prior model for the pressure-volume relationship. 
That is, while the constants in the model for the relationship may be 
determined from the data, the form of the model must not depend on the 
data to be used in fitting the relationship. In addition, the less 


accurately the hypothesized model describes this relationship, the less 


a 


valid will be the resulting uncertainty statements. That is, inaccuracies 
in the hypothesized model will lead to systematic differences between true 
and fitted curves. In order to obtain valid uncertainty statements, 
bounds for such differences must be determined, and added to the statisti- 
cal uncertainty as systematic error limits. 

The interiors of large processing tanks do not generally conform to 
idealized geometrical shapes, such as cylinders. Often, however, the tank 
can be considered to be composed only of segments for which an idealized 
model is a good Te cresentaedion. In this paper it is assumed that the 
tank is composed of a nee number, k+l, of distinct and known regions 
where the idealized relationship between the two variables pressure, p, 


and volume, v, is given by 


p = £(v) 
= g1(v) Equ<uw < 6] 
= go(v) Ey <v< & 
= 21.41 (v) GroS vy S Epa 


In addition, continuity of the relationship at the interior "knots" &;, 
i=1,...,k is required. 

In all that follows, eae and height refer to the portion of the 
tank above the bottom of the diptube used to measure pressure. The 
portion of the tank below that point is known as the heel, and is not 


treated in this paper. The pressure measured is the difference in 


pressure between the bottom of the diptube and a reference point at the 
top of the tank. 
The pressure-volume relationship can be ascertained from the follow- 

ing two equations. At height h the volume in the container is 

h 

v = f A(x) dx 

0 
where A(x) is the cross-sectional area at height x. Also, when the liquid 
height is h, p = hog, where p is the density of the homogeneous liquid 
and g is the acceleration due to gravity. Using these two equations one 


easily obtains 


p/pg 
v=f A(x) dx. 
0 


ov 1 

Thus = = --A(p/pg) and hence in areas of the tank where A(x) is 
P pg 

constant the volume-pressure relationship is a straight line. 


If A(x) is constant* in region i, as it obviously is for at least 


some regions of the tank shown in Figure 1, then 


We assume that each pressure measurement has a random error 


associated with it, and that these errors are independent and 


* If A(x) is not constant on an interval, then p is not a linear 
function of v on that interval. The program under discussion uses the 
B-spline basis when higher-order polynomial splines are required, 
because as pointed out in reference [4], the use of simpler representa- 
tions of polynomial splines may lead to numerical instability. 
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normally-distributed, with mean zero and constant variance rae (Recall 
that volume is assumed to be measured with no significant error.) 
Because of these errors, only estimates of the coefficients yy and 8j 
are obtained during the calibration process (experiment). These 
coefficient estimates are then used during plant operation to obtain 
estimates of the volume in the tank, utilizing the inverse of the 
relationship (1). 

Determination of uncertainty limits on these estimates is not 
erivial! "there are two sources of random error: estimation of the 
coeffit¢ients and 8 in the calibration experiment, and measurement of p 
during ‘operational’ use of the tank. The familiar linear regression 
model’ has properties, such as the nonexistence of means and variances of 
reciprocals, which make the analysis difficult. Special justification 
involving asymptotic (large sample-size) behavior is thus required in 
order to use a propagation-of-error approach to obtain appropriate 
approximate uncertainty limits on the estimated volumes. Furthermore, 
unless the p-v relationship is linear, normally distributed errors in 
the p-measurements during operation produce non-normal errors in the 
resulting estimates of v. The usual propagation-of-error technique does 
not take into’ account the differing characteristics of these errors. 
The new technique presented in this paper, in contrast to the 
propagation-of-error approach just mentioned, does allow a correct 
accounting for both. 

The calibration chart (i.e., the table of uncertainty limits) is 
produced after choosing two probabilities, a and 6. An exact statement 


giving the interpretation of these probabilities may be found in Scheffe 


[1] and in Scheffé, Rosenblatt, and Spiegelman [2]. However, an 
expanded, more heuristic explanation is given here. First, we require 
bounds for the calibration curve which will contain the entire curve 
with a prechosen probability 1-6. (Thus, 6 can be thought of as 
describing the uncertainty level to be associated with the outcome of 
the initial calibration experiment.) These bounds guarantee, with 
probability 1-6, that for any and every future volume v within the range 
of calibration of the tank, the v-interval (see Figure 2) that would be 
obtained by projection of the value f(v) through the curves to v and v 
would contain v. The second probability level to be chosen is a. (Here 
a can be thought of as describing the uncertainty level to be attributed 
to errors in future individual pressure measurements.) If o were known, 
we could state that the true pressure f(v) at the unknown volume v lies 
within the l-a confidence interval (p-zj-g/20, p + Zj-~q/29) with 
probability l-a, where p is the observed pressure and 2}j-,/2 is the 
two-sided l-a value for a normal distribution. The Scheffe procedure 
expands this interval appropriately, to account for the facts that o is 
estimated and that this estimate is used for the 1-6 bound on the curve 
and for bounds on many different f(v). It then combines the l-a 
confidence interval for f(v) with the 1-6 bounds on the calibration 
curve to produce calibration intervals I(p) for v. Construction of the 
calibration intervals is shown schematically in Figure 3. A set of 
intervals I(py) for py in the range of values obtained during the 
calibration experiment is called a calibration chart (see Figure 4). 

In the discussion of the example presented in the next section, more 


detail on the nature of the steps that make up a calibration run will be 


found. 


3. An Example. 

This example relates to a processing tank, roughly circular in cross 
section, une olin internal structure consisting of cooling coils, 
Pees praken, ete. [5]. This tank is pictured in Figure 1. The data 
from ealiiration runs on this tank have graciously been made available 
by the author of reference [5]. 

There were five calibration runs for which the data were useful for 
this Euivata: One run was done in the canyon where the tank is to be 
used. The other four, done in a mock-up area, used smaller tubing in 
the Preasnve-ncaguring system. This smaller tubing was known to cause 
Br onarie Aneees iene in the pressure measurement, which were expected 
to be Hneerly velated to pressure for each run. Since the tubing in 
the canyon was sufficiently large to render the pressure drop 
Paterieicant the vetenarie error was estimated for each of the other 
four runs, eau a correction made by applying a linear transformation to 
the measured pressure. It should be noted that these corrections, made 
to four of the five runs, effectively decrease the degrees of freedom 
for the error sum of squares by eight (two correction parameters times 
four runs corrected). 

The calibration program was applied to these data, as were various 
other ee eaatatien available ~ the large central computer at NBS. The 
results ea. ee be presented and their use described. 

In the version of the program described here, the knot locations are 
input by the apatineie. It is presumed that the knot locations can be 
adequately prescribed from the blueprints and other knowledge about the 


tank. (A refinement which allows the automatic determination of the 


number and location of knots is being investigated.) The program 
displays the given knot locations and other input data, as shown in 
Figure 5. Next come the results of the fitting operation, as shown in 
part in Figure 6. Note that the fit here is a fit of observed pressure 
(y) as a function of the accurately-dispensed volume (x); it is pressure 
which is subject to errors of observation, and volume which is 
essentially known. The residual standard deviation, an estimate of the 
standard deviation of the pressure measurements, is derived from the 
residuals or deviations of the measured pressures from the fitted curve. 
In this case, its value is 1.49 pascals. This value, the corresponding 
degrees of freedom, and the coefficients (which are in general not 
immediately interpretable, since they refer to the so-called B-splines, 
a representation chosen for computational stability), are part of these 
results. The calibration intervals for estimation of volume from 
measured pressure are printed next (see Figure 4). An ordinary 
polynomial representation and a residual plot are also printed as shown 
in Figures 7 and 8. 

It will be instructive to examine the printout and discuss the 
approach in more detail, and this will now be done. 

As can be seen from Figure 5, the program duplicates the end pisces 
This is just a simple way to define the B-spline basis functions which 
are used to perform the fit, and need not concern the analyst. The 
input values for knot locations, degree of fit, and other miscellaneous 


parameters are printed out for verification. 


aw ee EE 


At this point, the program does a linear least squares fit of the 
specified model to the (v,p) data, and prints out a reasonably standard 
summary (Figure 6). The column labels are self-explanatory. At the 
bottom of this summary are found the residual standard deviation and its 
associated degrees of freedom, and the estimated coefficients with the 
corresponding estimated standard deviations. 

The program next computes some intermediate results which generally 
are of no interest to the analyst, and therefore are only printed out if 
requested. These are confidence intervals for p, at 300 evenly-spaced 
points covering the range of v between the extreme knots. Input values 
of a, 6, and c are used in this procedure, so these values are printed. 

The calibration chart comes next, giving the predicted value of v 
and the corresponding lower and upper limits for each of the specified 
set of p-values (see Figure 4). It is obtained by inverse interpolation 
from the confidence intervals for p discussed in the preceding 
paragraph. Usually, the extreme values of p will be at least partially 
outside the range of at least one of the curves. When this happens, the 
intervals should extend either to zero volume or to full volume. This 
is indicated by "<" and ">" respectively on the printout. 

Since the coefficients actually fitted are the B-spline coefficients, 
the program converts the B-spline representation to a simple polynomial 
representation. The printout shows the endpoints and the coefficients of 
the fitted polynomial for each of the specified intervals (see Fig. 7). 

Finally, the residuals from the fitted model are plotted in order of 
increasing volume to allow a visual check of the adequacy of the 


chosen model [6] (see Fig. 8). 


At NBS, with the aid of the central computer and the OMNITAB system 
[7], a number of other things were tried which strengthen the conviction 
that this program does indeed work well. These will now be discussed. 

Various subsets of the data were fitted to the same model. No incon- 
sistencies were found. 

The sensitivity to position and presence of the different knots was 
checked. The results were rather sensitive to the knot locations, which 
implies that good estimates of the locations are required for good fits. 
It should also be noted that where a knot bounds a short segment, the 
removal of that knot might make very little difference in any global 
measure of fit quality, unless there are many data points in that short 
stretch. Nevertheless, the systematic error introduced by deleting that 
knot can be a consistent source of inventory losses or gains, apparent or 
real. Thus it is important to include all real segments in the model to be 
fitted. 

A separate program was written to perform linear spline fitting, while 
the main package was being put together. The answers did not differ 
between the two programs, providing a partial check that no programming 
errors were committed. 

Smooth higher-order spline fits were tried (quadratic and pee 
There was no improvement in fit. The linear spline appears to provide 
an adequate representation of the pressure/volume ee a 

Probability plots were done in various ways, looking for possible 


troubles with the data or the method. Nothing suspicious was found. 
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4. User's Manual. 

This fixed-knot spline package for calibration consists of a "main" 
subroutine SPLEEN and 29 additional subroutines. The manner in which 
they interact is diagrammed in Figure 9. All programs are written in 
FORTRAN and have been checked for portability by the Bell Laboratories 
PFORT verifier [8]. It was decided that SPLEEN should be a subroutine 
rather than a main program so that the user could enter the parameter 
values in the way most convenient for him. The user then must write a 
main program which sets up the required dimensioned variables and 
assigns values to cha) Geaeares parameters (those with asterisks in the 
list which follows). These parameters are passed to subroutine SPLEEN 
via the statement 

CALL SPLEEN(H,X,Y,W,R1,R2,RES,N,NX,NKX,T, BCOEF ,XX1,Q,DIAG,K, 
KX, YY,NY,NYX,MD,SCRTCH,JX,AL,DL,C, IP) 


where 


* H(80) = Up to 80 characters in 80Al format identifying the 
‘data 


* X(NX) = Vector (length N) of X-values where observations were 
made (independent variable) 


* Y(NX) = Vector (length N) of observations 
* W(NX) = Vector (length N) of weights for observations 
R1(NKX) = Vector (length N+K) for scratch area 
R2(NKX) = Vector (length N+K) for scratch area 
RES(NKX) = Vector (length N+K) of residuals from spline fit 
* N = Number of observations 
* NX = Dimension (AN) of vectors X,Y,W 
* NKX = Dimension (NK) of vectors R1,R2,RES 


* T(KX) = Vector (length K+2*MD) of knot locations 


at 


BCOEF (KX) 


XXI (KX,KX) 


Q(JX, KX) 
DIAG(KX) 


* K 


* KX 


* YY (NYX) 


* NY 


* NYX 


* MD 


SCRTCH(JX,JX) 


* JX 


* AL 


Vector (length K+MD-1) of B-spline coefficients 


Variance-covariance matrix (size [K+MD-1]x[K+MD-1] ) 
of B-spline coefficients 


Matrix (size [MD+1]x[K+MD-1]) for scratch area 
Vector (length K+MD-1) for scratch area 


Number of knots specified by user (later increased to 
K+2*MD by program) 


Dimension (?K+2*MD) of vectors T,BCOEF,DIAG and 
matrices XXI and Q (column only for Q) 


Vector (length NY) of Y-values for which predicted 


X-values (with confidence intervals) are to be 
computed 


Number of Y-values for which predicted X-values are 
to be computed 


Dimension (2NY) of vector YY 


Degree of spline (<19); for example, l=linear, 
2=quadratic, 3=cubic) 


Matrix (size [MD+l1]x[MD+1]) for scratch area 


Dimension of square matrix SCRTICH and row dimension 
of matrix Q = 20 


Alpha level of significance 
Delta level of significance 


Constant in the interval (0.85,1.25) associated with 
Scheffe's calibration technique (see Appendix 1 for a 
discussion of this constant) 


OQ For abbreviated printout 

1 For full printout (residuals, PP representation) 

2 For full printout plus Y-confidence intervals for 
300 evenly spaced X-values over knot span 


Variables which appear with an asterisk (*) require input values 


from the main program. The subscripts on vectors and matrices 


indicate the dimensions which must be assigned in the main program. 
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Variable names which begin with the letters I1,J,K,L,M, or N are of the 
INTEGER type. The remaining variable names are of the REAL single 
precision type. 

The print parameter IP gives the user a certain amount of control 
over the amount of information to be printed out. Normally the most 
suitable value is IP=l. A value of IP=0 suppresses the printout of the 
weights, independent variable, observations, predicted values, and 
residuals. This option may save quite a bit of paper in case there are 
several hundred observations, but it deprives the user of the chance to 
visually examine the residuals. A value of IP=2 causes a listing of 
certain intermediate vectors which are somewhat lengthy and would not 
normally be of use to the user. 

In the interest of minimizing the number of variables needed in 
the CALL statement, not all of the printed information can be recovered 
through the passed parameters. Furthermore, three of the variables (X, 
K, and T) return values different from their input values. 

The data points (X;, Y;, Wy) may be input in arbitrary order, as may 
the knot locations Ty and the vector of YY; specifying the y-values on 
the calibration chart. 

There are two subroutines which check for conediecency among the 
input parameters. Each inconsistency causes a diagnostic message to be 
printed. If one or more inconsistencies is detected then the program 
execution is terminated. Observations outside the knot span are flagged 
and weighted zero. The number of observations is then reduced by one 
for each flagged point and a diagnostic is printed. This is not a fatal 


error unless it reduces the number of degrees of freedom to zero or less. 


as) 


Although this package can handle splines of any degree up to 19 it 
was primarily intended for splines of lower degree, i.e., linear, 
quadratic, or cubic. Test runs on sets of both real and artificial data 
have given valid results up 88 about degree 9. Beyond that the 
limitation of single precision arithmetic on the 36-bit NBS central 
computer begins to cause roundoff errors that invalidate the results. The 
user should exercise caution when fitting the higher degree splines. 

If the user wants to change some of the continuity conditions at a 
given knot he may do so by duplicating that knot in the knot vector 
which is passed to subroutine SPLEEN. If a knot appears M times in the 
fitting of a spline of degree N then the functional value and the first 
N-M derivatives of the function will be continuous at that knot. If 
M = N+l, neither the function nor its derivatives are required to be 
continuous. 

The package may be applied to both monotone increasing and 


decreasing calibration curves. 


5. Summary and Discussion. 

An approach to calibration curves and their uncertainty bands has 
been presented, complete with a FORTRAN program to perform the required 
calculations. An example involving a large process tank has been used to 
illustrate the approach and the program. The results include not only 
the curve for estimating volume from measured pressure, but also valid 
uncertainty limits for repeated applications of the calibration curve 
obtained. 

The interval estimates of volume comprise two parts: a long-term 


component which changes only at recalibration, and a random component. 
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The contribution due to the long-term component may be estimated in 
large scale calibration experiments by the volume interval estimate 
obtained set aie Similarly, the contribution due to the random 
component may be estimated by the volume interval estimate obtained when 
6=1. When the calibration experiment is of a more modest size involving 
less than 100 ae atnee the above component estimates may not be 
realistic. However a more comprehensive treatment for combining 
interval estimates (and hence their components) obtained from a calibra- 
tion curve is under development by C. Spiegelman and K. Eberhardt [9]. 

The results of a calibration will be used repeatedly, usually 
without any further opportunity to verify their correctness, until the 
next pal tiracion, Berek oee it is important that the measurement system 
be under control. In the werk reported here, the run-to-run differences 
observed in the mock-up area were due to a known source (the small 
diameter of the tubing), and could be corrected. If any anomalous 
behavior is observed which cannot be satisfactorily explained, then of 
course the entire statistical analysis must be approached with caution. 

Little has been written about the design of calibration experiments 
- i.e., the Bal aebion of volumes at which pressure is to be measured, 
the number of measurements to be taken at each volume, and the arrange- 
ment of these measurement points into a sequence of runs. One solution 
to this question has Hea achieved by Spiegelman and Studden, and will 
be published in the NBS Journal of Research [10]. In general, later runs 
will concentrate on west ese of a tank, but it is good practice 
to ensure that at ie tus runs cover each section, and that several 


runs cover the whole tank. If this precaution is not taken, there might 
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be very little cross-validation between runs. 
Certain caveats ought to be mentioned here. The program under 


discussion assumes that the knot locations are known, and that the model 


is correct. Consequences of failure of these assumptions could be 
severe. With respect to the knot locations, careful inspection of the 
residual plots will sometimes indicate discrepancies. These may be 
small; however, it is important to realize that such regions represent 
systematic deviations, and could be used (at least in theory) to cover 
the diversion of material. An approach to the problem involving unknown 
knot locations is being pursued at this writing. 

Unlike a simple linear regression, where the inclusion of a 
superfluous higher-order term generally causes no major trouble (the 
fitted coefficient turns out insignificantly small, and the residual mean 
square increases minimally), choosing a higher-order model when fitting 
smooth splines can result in a very much worse fit. This is because of 
the smoothing restrictions, which greatly limit the freedom of the 
fitting procedure. (Imagine that the true model consists of two straight 
lines meeting at a point. If one chooses to fit a quadratic spline, then 
one is insisting on having two quadratic curves which meet at the proper 
x-value, and which have the same slope at that point. Thus the slope at 
that point is probably going to be some value between the two 
straight-line slopes, and the fit cannot be accurate.) One way around 
this difficulty is to fit piecewise polynomials (i.e., do not require 
smoothness), and investigate the appropriate degree from these fits. 
However, it is much better to know the situation well enough to choose 


the correct model from physical considerations. 
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A technique that the authors have found useful is to run the program 
described here with the degree of the fit set equal to zero. The result 
is to fit a step function to the Bees and to produce a plot of the 
residuals from that fit; for the Example of this paper, that plot is 
repeeaieen as Figure 10. It one SE seen that the residuals look as 
linear as a printer plot eericek. Therefore, a first-degree (linear) 


spline fit is the proper choice. If in some segment the relationship 


were not linear, this plot should show ite The plot also gives some idea 


of the spread of points across the intervals, though of course near~- 


duplicate points will plot as one because it is a discrete printer plot. 
Note that the continuity restrictions can be relaxed when using the 
program under consideration, by simply duplicating the knots. See 


Section 4 for details. 


wy 


Appendix 1. 


As stated in Section 4, a constant c must be input by the user. In 
order to obtain this constant from tables 1 and 2 in Scheffe (1973) for 
1<ps<l0 the oe must have calculated the standard deviation, SD, for p(v) 
in the complete region of calibration. The smallest and largest values 
the sp(p(v))/o over the complete calibration region are used as input for 
the Scheffé tables. For k>10 Scheffe gives a mathematical algorithm for 
finding c, and states that for very large (asymptotic) values of n-k, c=l. 
(Here n is the number of observations, and k is the number of B-spline 
coefficients.) If the reader does not wish to do a Scheffe table 1 or 
table 2 lookup, the following table gives approximate and generally larger 
values for Gife omastent. 

Approximate c values 


for 1<k<l0 


n-k 60-119 120-149 150 + 
c 1.10 1.05 1.00 
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Appendix 2. 
Program Listing. 
The subroutines which make up the spline-fitting package follow, in 


alphabetical order. 
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* COMPUTATICN OF CALIBRATION INTERYALS BY INVERSE INTERPOLATION * 
LOWER LIMIT PREDICTED UPPER LIMIT 
I YC I) FOR X x FOR X 
1 3713.009 1.705250 1.705890 1.769779 
2 3823 . 600 1.749691 1.733002 1.756223 
3 3933. 000 1.797422 1.800113 1.892745 
4 4943 .000 1.844952 1.647225 1.849473 
5 4153 .000 1.892690 1.894297 1.696511 
6 4263 .960 1.939062 1.941202 1.943401 
ré 4373 .000 1.985893 1.988067 1.999239 
& 4483 .690 2.632763 2.034931 2.037977 
9 4393 .000 2.979673 2.081796 2.083916 
10 4703 .000 2. 126563 2. 1286690 2.130755 
164 21643.60 9.704492 9.796814 9.709139 
165 21753.00 9.755669 9.738022 9.760378 
166 21863 .00 9.800846 9.899230 9.611617 
167 21973.90 9.858923 9.860439 9.862257 
168 22983 .00 9.999199 9.911647 9.914998 
169 22193.00 9.960375 9.962855 9.965333 
170 22393 .60 1@.01155 10.01466 10.61656 
171 22413.60 10. 06267 10.06506 10.@6742 
172 22523. 80 10.11361 10.11578 10.11793 
173 22633 .90 10. 16435 10. 16651 10. 16866 
6 
234 29343 .60 13.28631 13.28856 13.29082 
235 294353 .00 13.33740 13.33971 13.34263 
236 29363 .90 13.38549 13.39@87 13.39325 
230 29673 .00 13.43957 13.44262 13. 44447 
236 29783 .99 13.49066 13.49317 13.49570 
239 29393 .60 13.54174 13.34433 13.54692 
240 39003 .09 13.59251 13.59548 13.59315 
241 39113.00 13.64334 - 8996699 > 13.64334 
242, 39223. 09 13.64334 - 6809099 > 13.64334 
*xxk AT LEAST ONE Y-VALUE IS OUTSIDE THE RAGE 7: 
RK GF AT LEAST ONE CALIBRATION CURVE se oKIR 
S DENOTES THE VALUE CF THE SMALLEST KNOT 
L, DENOTES THE VALUE OF TH LARGEST KNOT 


* DENOTES VALUES OUTSIDE THE RANGE OF THE CALIBRATION DATA — NO 
VALID PREDICTION IS AVAILABLE 


< INDICATES THAT NO VALID LOWER CALIBRATION LIMIT 
GREATER THAN THE MINIMUM POSSIBLE X-VALUE IS AVAILABLE. 


> INDICATES THAT NO VALID UPPER CALIBRATION LIMIT 
SMALLER THAN THE MAXIMUM POSSIBLE X-VALUE IS AVAILABLE. 


; 
Figure 4. Calibration chart. Y is pressure in pascals; X is volume in™M. 


95 


SERRA NRE RR I RR \ 
x FIXED-KNOT SPLINE PACKACE FOR CALIBRATION = 


HANAN ANE RARER RE REAL NE EIR aq 
REAL DATA FROM A TANK CALIBRATION (DATA IN FILE CPReNS.95) SI 
a 

----- FULL PRINTOUT ----- . 2 

) | A 

<<<<< EACH END KNOT DUPLICATED 1 TIMES >>>>> fl 


ee ew ae a es ee SS eS 


KNOTSC TI) 


I 

1 1.70525 | 
2 1.70525 

3 1.89466 
4 4.54674 
5 5.49482 
6 5.87331 q 
ra 6.44156 2 7 
re 7.18990 
9 10.0425 
10 10.2320 
11 12.5044 a 
12 13.6483 | 
13 13.6483. ; 4 


re es sw Se me ee ee es ee ae er ee ee Oe ee eee 


DEGREE OF SPLINE = 1 
NUMSER OF OBSERVATIONS = 172 | 
NUMBER OF ZERO WEIGHTS = 1) 
NU?HSER OF NON-ZERO WEIGHTS = 172 
NUHBER OF KNOTS = 13 
NUMBER OF B-SPLINES = 11 


NUIBER OF Y VALUES 
FOR WHICH X COMFIDENCE = 
INTERVAL IS TO BS COMPUTED 


{<<< 11 Be-SPLINE COEFFICIENTS COMPUTED >>>>> ! 


il 
rN) 
+>» 
to 
PSS © rr rans ee regmarc sree comn 


Figure 5. Preliminaries. 
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= 


HSSatagaA Wh = 


—_ 


162 
163 
164 
165 
166 
167 
168 
169 
170 
171 
172 


RESIDUAL STD DEV 


1.9000 
1.0000 
1.9000 
1.0000 
1.0006 
1.0006 
1.9000 
1.0000 
1.0000 


- 9000 
- 9000 
- 0600 
- 9000 
- 0000 
1.92000 
1.9900 
1.9009 
1.69000 
1.9900 
1.90000 


et ete 


1.48848 


XC 1) 


1.795250 
1.705270 


. 1.894010 


1.894140 
1.894580 
1.894650 
1.894660 
2.084060 
2.084110 
2.272720 
2.273020 


13.07486 . 


13.26022 
13.26158 
13.26232 
13.45198 
13.45385 
13.639 14 
13.640352 
13.64133 
13.64143 
13. 64333 


161 


OBSERVED 
YC 1) 


3711.640 
3711.420 
4152.240 
4152.130 
4151.050 
4151.910 


4153 .460. 


4599 .320 
4600 .990 
5044. 100 
5042.730 


28884. 68 


29282 .92 
29285.71 
29284.50 
29695. 29 
29699. 14 
30096.50 
30099 .62 
30098. 26 
30102.36 
30105.91 


RESIDUAL D.F. 


i a se es a ee ee ee ee ee 


= 


HROvOAninadhoane 


—a 


K<<K< 


B-SPLINE 
COEF 


3711.5053 
4153 .7559 
10378. 705 
12579 .214 
13403. 155 
14630.321 
16236 . 443 
22364. 0a4 
22764. O45 
27656 . B46 
30105 .922 


STD DEV 


STD. DEV. 


1.0325651 
- 41250259 
41074148 
- 37466792 
- 4306913 
- 62254084 
-41161732 
- 44698178 
- 43252287 
- 39959609 
- 32032419 


PREDICTED 
YC) 


3711.505 
3711.552 
4152.238 
4152.542 
4153 .569 
4153.732 
4153.756 
4598.315 
4598. 432 
59041. 136 
5041.841 


28863 . 46 
29282. 06 
29284.98 
29286 .57 
29694. 42 
29698. 44 
30096 .89 
30999. 86 
30101.60 
301061.61 
30105.90 


<<<<< GRID OF 300 EVENLY SPACED X VALUES CREATED >>>>> 


Figure 6. Regression Bit. 
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RES IDUAL( I) 


- 13455 
-, 13223 
- 18311-92 
-.41168 
-2.5190 
-1.6225 
-—.35577 
1.0054 
1.6581 
2.9636 
- 88947 


1.2214 
- 86230 
- C2754 
-2.90735 
- 67061 
- 69897 
—.38940 
—.23766 
-3.3391 
- 346 14 
- 10254-61 


OF 360 PREDICTED Y VALUES COMPUTED >>>>> 


STD DEV OF 
PREDICTED YC I) 


1.052565 
1.052453 
-4110981 
-4113765 
- 4123280 
- 4124867 
-4125026 
3729 185 
- 3729083 
- 3359 134 
- 3358569 


- 2921614 
- 3425575 
- 3430482 
- 3433159 
- 423576 1 
- 4244610 
-3160907 
- 9186241 
-9192548 
- 3193080 
- 5203189 


SLLOL_ILLI_L LSD SS ND SS SS SS SS SS SD Se SSS SD a Oe Se Ow Oe Oe Os Soe ye ee oe ees oe ese wo ee ee we 


esse INTERVAL..... COEFFICIENTS OF (X=NC1))+*P 
I XC) XC I+1) P= 6 i 
1 1.7052 1.8947 3711.5 2334.9 
2 1.8947 4.3467 4152.8 2547.2 
3 4.5467 . 5.4948 10379. 2321.0 
4 5.4948 5.8738 12579. 2174.0 
5 5.8738 6.4416 13403. 2161.5 
6 6.4416 7.1899 14639. 2146.2 
4 7.1899 10.043 16230. 2146.1 
8 10.043 10.232 22364. 2168.6 
9 10.282 12.504 22775. 2148.3 
10 12.504 13.643 276357. 2150.4 
. Figure 7. Ordinary polynomial 


representation of the Fitted spline. 
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Figure 8. 
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Figure 9. Diagram of subroutine interactions. 
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